Study of the Decays B° -» £>(*)+£>(*)- 



E. Lipeles, M. Schmidtler, A. Shapiro, W. M. Sun, A. J. Weinstein, and F. Wiirthwein* 
California Institute of Technology, Pasadena, California 91125 

D. E. Jaffe, G. Masek, H. P. Paar, E. M. Potter, S. Prell, and V. Sharma 
University of California, San Diego, La Jolla, California 92093 

D. M. Asner, A. Eppich, T. S. Hill, R. J. Morrison, and H. N. Nelson 

University of California, Santa Barbara, California 93106 

R. A. Briere 

Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 

B. H. Behrens, W. T. Ford, A. Gritsan, J. Roy, and J. G. Smith 

University of Colorado, Boulder, Colorado 80309-0390 

J. P. Alexander, R. Baker, C. Bebek, B. E. Berger, K. Berkelman, F. Blanc, V. Boisvert, 
D. G. Cassel, M. Dickson, P. S. Drell, K. M. Ecklund, R. Ehrlich, A. D. Foland, 
P. Gaidarev, L. Gibbons, B. Gittelman, S. W. Gray, D. L. Hartill, B. K. Heltsley, 
P. I. Hopman, C. D. Jones, D. L. Kreinick, M. Lohner, A. Magerkurth, T. O. Meyer, 
N. B. Mistry, E. Nordberg, J. R. Patterson, D. Peterson, D. Riley, J. G. Thayer, 
P. G. Thies, B. Valant-Spaight, and A. Warburton 
Cornell University, Ithaca, New York 14853 

P. Avery, C. Prescott, A. I. Rubiera, J. Yelton, and J. Zheng 
University of Florida, Gainesville, Florida 32611 

G. Brandenburg, A. Ershov, Y. S. Gao, D. Y.-J. Kim, and R. Wilson 

Harvard University, Cambridge, Massachusetts 02138 

T. E. Browder, Y. Li, J. L. Rodriguez, and H. Yamamoto 
University of Hawaii at Manoa, Honolulu, Hawaii 96822 



T. Bergfeld, B. I. Eisenstein, J. Ernst, G. E. Gladding, G. D. Gollin, R. M. Hans, 
E. Johnson, I. Karliner, M. A. Marsh, M. Palmer, C. Plager, C. Sedlack, M. Selen, 

J. J. Thaler, and J. Williams 



University of Illinois, Urbana- Champaign, Illinois 61801 



K. W. Edwards 
Carleton University, Ottawa, Ontario, Canada K1S 5B6 
and the Institute of Particle Physics, Canada 

R. Janicek and P. M. Patel 

McGill University, Montreal, Quebec, Canada H3 A 2T8 
and the Institute of Particle Physics, Canada 

A. J. SadofT 
Ithaca College, Ithaca, New York 14850 

R. Ammar, A. Bean, D. Besson, R. Davis, N. Kwak, and X. Zhao 
University of Kansas, Lawrence, Kansas 66045 

S. Anderson, V. V. Frolov, Y. Kubota, S. J. Lee, R. Mahapatra, J. J. O'Neill, R. Poling, 

T. Riehle, A. Smith, and J. Urheim 

University of Minnesota, Minneapolis, Minnesota 55455 

S. Ahmed, M. S. Alam, S. B. Athar, L. Jian, L. Ling, A. H. Mahmood, 1 " M. Saleem, 

S. Timm, and F. Wappler 

State University of New York at Albany, Albany, New York 12222 

A. Anastassov, J. E. Duboscq, K. K. Gan, C. Gwon, T. Hart, K. Honscheid, D. Hufnagel, 
H. Kagan, R. Kass, T. K. Pedlar, H. Schwarthoff, J. B. Thayer, E. von Toerne, and 

M. M. Zoeller 

Ohio State University, Columbus, Ohio 43210 

S. J. Ri chichi, H. Severini, P. Skubic, and A. Undrus 
University of Oklahoma, Norman, Oklahoma 73019 

S. Chen, J. Fast, J. W. Hinson, J. Lee, N. Menon, D. H. Miller, E. I. Shibata, 

I. P. J. Shipsey, and V. Pavlunin 

Purdue University, West Lafayette, Indiana 47907 

D. Cronin-Hennessy, Y. Kwon,* A.L. Lyon, and E. H. Thorndike 



2 



University of Rochester, Rochester, New York 14-627 

C. P. Jessop, H. Marsiske, M. L. Perl, V. Savinov, D. Ugolini, and X. Zhou 

Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309 

T. E. Coan, V. Fadeyev, Y. Maravin, I. Narsky, R. Stroynowski, J. Ye, and T. Wlodek 

Southern Methodist University, Dallas, Texas 75275 

M. Artuso, R. Ayad, C. Boulahouache, K. Bukin, E. Dambasuren, S. Karamov, 
G. Majumder, G. C. Moneti, R. Mountain, S. Schuh, T. Skwarnicki, S. Stone, 
G. Viehhauser, J.C. Wang, A. Wolf, and J. Wu 

Syracuse University, Syracuse, New York 13244 

S. Kopp 

University of Texas, Austin, TX 78712 

S. E. Csorna, I. Danko, K. W. McLean, Sz. Marka, and Z. Xu 
Vanderbilt University, Nashville, Tennessee 37235 

R. Godang, K. Kinoshita, § I. C. Lai, and S. Schrenk 

Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

G. Bonvicini, D. Cinabro, S. McGee, L. P. Perera, and G. J. Zhou 

Wayne State University, Detroit, Michigan 48202 

(CLEO Collaboration) 
(February 7, 2008) 



*Permanent address: Massachusetts Institute of Technology, Cambridge, MA 02139. 
^Permanent address: University of Texas - Pan American, Edinburg, TX 78539. 
^Permanent address: Yonsei University, Seoul 120-749, Korea. 
^Permanent address: University of Cincinnati, Cincinnati, OH 45221 



3 



Abstract 



The decays B° -» D*+D*-, B° -» D^D^ and £° -► D+D" are stud- 
ied in 9.7 x 10 6 T(4S) -> SB decays accumulated with the CLEO detec- 
tor. We determine £(5° -> D* + D*-) = (9.9±f § [stat.] ± 1.2 [syst.]) x 10" 4 
and limit -» D^D^) < 6.3xl(T 4 and B(B° -» D+D") < 9.4X10" 4 at 

90% confidence level (CL). We also perform the first angular analysis of the 
B° — > D* + D*~ decay and determine that the CP-even fraction of the final 
state is greater than 0.11 at 90% CL. Future measurements of the time de- 
pendence of these decays may be useful for the investigation of CP violation 
in neutral B meson decays. 
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Typeset using REVTpjX 



I. INTRODUCTION 



The first observation of CP violation outside the neutral kaon system |T]|2| may well be 
a non-zero difference in the rates of P° — > J/ipK§ and P° — > J/ipK® decays 0. Such a 
measurement would be an important test of the Standard Model mechanism for CP violation 
as described by the CKM quark- mixing matrix |4j . In the Standard Model, the CKM matrix 
is unitary; for three quark generations, this property can be represented as a triangle in the 
complex plane with internal angles a, (3 and 7 ||. Asymmetries in the rate of neutral B 
meson decays to CP eigenstates that occur via the Cabibbo-favored b — > cW + ; W + — > cs 
(eg., B° — > J/ipKg) process are expected to be proportional to sin 2/3. In contrast to the 
decay B° — > J/ipK^, for the Cabibbo-suppressed processes [] B° -f £>(*)+£>(*)-, the weak 
phase difference between the tree (b — ► ccd) and penguin (b — ► dec) amplitudes may be 
appreciable [0,0. In the absence of a strong interaction phase difference between the tree 
and penguin B° — ► p(*)+p(*)~ amplitudes, the magnitude of the asymmetry would be also 
proportional to sin 2(3. The decay rate asymmetry of P° — > D* + D*~ decays would also be 
proportional to sin 2(3 but may suffer from dilution due to the P-wave (CP-odd) component 
of the D* + D*~ final state The relative CP-even and CP-odd components of the 

B° — > D* + D*~ decay can be determined by an angular analysis jlO| that removes any such 
dilution. 

Measurements of rate asymmetries in the decays B° — > pM+p(*)~ may provide a 
means to resolve the four-fold ambiguity in (3 inherent in a measurement of sin 2(3 from 
P° — ► J/tpKg decays ||8| JTT|JT2]| . Comparison of the measured asymmetries in B° — > J/ipK^ 
and B° — > D + D~ decays may allow partial resolution of the ambiguity in the determination 
of (3 if the sign of the ratio of the tree and penguin amplitudes of B° — > D + D~ decays can be 
ascertained [|13|]. B° and B° mesons decay to the same D* + D~ final state with amplitudes 
of comparable magnitude and significant interference between them is possible |Hyi5(| . As 
for B° — > D + D~, the asymmetry between the rates of B° — > D* + D~ and B° — > D*~D + is 
directly proportional to sin 2(3 in the absence of strong phase differences. In the presence 
of a strong phase difference, the rate asymmetry would depend on both sin 2(3 and cos 2(3 
and, when combined with a sin 2(3 measurement from B° — > J/i/jK® decays, could aid in the 
resolution of ambiguities in the determination of (3. 

The decay B° — ► P* ± P =F would also provide a clean test of the factorization ansatz for 
decays into two charm mesons and provide a measurement of the ratio of D* + and D + decay 
constants and form factors P,|l4|. 

The expected branching fractions of the decays B° — ► p(*)+p(*) _ can be estimated from 
the measurement of the corresponding Cabibbo-favored processes B° — > pW+pM - || and 
the ratio of decay constants []8|,|IE|. The estimated P° — > D* + D*~ branching fraction is 
~ lOxlO -4 , consistent with the measurement of (6.2^9 [stat.] ± 1.0 [syst.]) x 10 -4 ||17|| ) and 
the estimates for P° — > D* ± D T and P° — > P + P _ are ~ 8xl0~ 4 and ~ 5xl0~ 4 , respectively. 
We present an update of the previous CLEO measurement of B(B° — > P* + P* _ ) ]1^| and 



improved upper limits on i3(P° — > P* ± P =F ) and i3(P° — > P + P ) |T3 based upon a sample 



1 B° dM+dM- denotes the decays B°^D+D-, B° -» D* + D~, B° -> D + D*~ and 
50 _^ D*+D*~. B° -> P*±PT denotes the sum of P° P*+P" and S° -» D + D*-. 
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of 9.7 x 10 6 BB pairs produced in e + e _ ^T(4S) decays accumulated with the CLEO detector 
at the Cornell Electron Storage Ring (CESR). We also present the first angular analysis of 
B° — > D* + D*~ decays and limit the CP-odd content of this reaction. The results presented 
here supersede the previous CLEO results |T7|Jl8| . 



II. THE CLEO DETECTOR 



The data were accumulated with two configurations of the CLEO detector dubbed 



CLEO II [[HJ and CLEO II. V [pOj] . In the first configuration, a 1.5T solenoidal magnetic 
field encloses three concentric cylindrical drift chambers that are nested within a cylindrical 
barrel of time-of-flight (TOF) scintillators and a CsI(Tl) calorimeter. The surrounding iron 
return yoke is instrumented with proportional wire chambers for muon identification. The 
large outer drift chamber provides up to 49 measurements of a charged particle's specific 
ionization (dE/dx) for particle species identification. In the CLEO II. V configuration, the 
innermost wire chamber was replaced by a three-layer, silicon vertex detector (SVX) capable 
of providing precision position information in both rcf) and z . The gas in the large outer 



drift chamber was also changed from argon-ethane to helium-propane, resulting in improved 
dE/dx and momentum resolution [22 . 



The Monte Carlo simulation of the CLEO detector response was based upon GEANT p3 



Simulated events for the CLEO II and CLEO II. V configurations were processed in the same 
manner as the data. 



III. CHARM MESON RECONSTRUCTION 

Observation of the relatively small rates expected for B° — > decays requires 

an aggressive program of charm meson reconstruction. The D° decay modes considered for 
reconstruction are K~rr + , K~7r + 7r°, K^tt+tt+tt^ , Ag7r + 7r~ and A"gV + 7r~7r ; the D + decay 
modes considered for reconstruction are K~7t + tt + , Ag7r + , Ag7r + 7r°, A"gV + 7r + 7r~, K~7t + 7t + 7t , 
K~K + tt + and K~ K + tt + tt . In order to limit background, the decay D + — > K~7t + 7t + 7t° 
is not considered for the reconstruction of the B° — > D* + D*~ mode and the decays 
D + — > K~K + tt + and D + — > K~ K + tt + tt° are not considered for the reconstruction of the 
B° — > D + D~ mode. The D* + decays to D°tt + and D + ir° are selected for the reconstruction 
of the B° -> D*+D*- and B° D*±DT modes, although the final state (D + tt°)(D-tt°) is 
overwhelmed by combinatorial background and is excluded from the B° — > D* + D*~ recon- 
struction. In the following, "D" refers to either D + or D° mesons, "7r s " refers to the slow 
pion daughter of the D* + decay and charge conjugation is implied unless explicitly stated 
otherwise. 

Charged kaon and pion daughters of D meson candidates must be compatible with an 
origin at the e + e~ interaction point. The dE/dx or TOF measurement of a charged track, 
when available, must be within 2.5 and 3.0 standard deviations (a) of expectations for K ± 
and 7r candidates, respectively. The K§ meson candidates are reconstructed in the tt~tt + 
decay mode and must be consistent with an origin at the e + e~ interaction point. At least one 
of the K§ daughter pions must be inconsistent with an origin at the e + e~ interaction point. 
Neutral pion candidates are formed from energy deposits in the calorimeter consistent with 
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electromagnetic showers unassociated with a charged track and with an energy exceeding 30 
MeV in the barrel (| cos#| < 0.71) and 50 MeV in the endcap region where 9 is the angle of 
the shower with respect to the z axis. A requirement on the ir° minimumn momentum of 100 
MeV/c is imposed for D daughter candidates and of 70 MeV/c for D* + daughter candidates. 
The charged and K$ daughters of all D meson candidates are required to originate from a 
common vertex. 



IV. B° MESON CANDIDATE SELECTION 

A number of observables are used to suppress backgrounds. In general, the require- 
ments on these are more stringent for the B° — ► D*^D T and B° — > D + D~ modes than 
for B° — > D* + D*~ because the combinatorial backgrounds are larger. In addition, while 
common selection criteria for all D* + and D decay modes of each B° candidate were sat- 
isfactory for the B° — > D* + D*~ mode, the B° — > £)* ± D =F and B° —> D + D~ modes require 
separate criteria for each B° — > D^ + D~ channel to reduce background. The selection 
criteria for each channel of the B° — > D and B° — > D + D~ modes were optimized 

using simulated signal and background events assuming B(B° — > D* ± D T ) = 8xl0~ 4 and 
B(B° -> D+D~) = 4.5 xlO" 4 , respectively. 

A. B° meson candidate energy and mass 

The observable AE = E(D^ + ) + E(D^~) - E heam exploits energy conservation for 
B°B° meson pairs produced in T(4S) decays and has a resolution a(AE) = 8 MeV after 
constraining the B° daughter candidates to the masses [0. The beam-constrained B 

mass is defined as M(B) 2 = E^ eaum — pg, where p# is the measured B° candidate momentum. 
The M(B) resolution of 2.5 MeV is dominated by the beam energy spread |24" |. Signal 



candidates are selected by requiring both AE and M(B) — M% to be within 2.5a of zero for 
the B° -> D* + D*- mode and within 2.0a of zero for the B° -> D* ± D T and B° -> D + D- 
modes, where Mg is the world-average B° mass 0. 

B. Candidate mass x 2 

The overall deviation of D* + and D candidates from the D* + and D meson masses is 
quantified by 

v 2 _ v / M-M- \ 2 ( AMj - AM° \ 2 



2 For example, for the B° — > D + D mode, there are a total of 30 possible channels for the five 
D + decay modes in each detector configuration. 
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where Mj is the measured D candidate mass, AMj is the mass difference between the D* + 
and D candidates, and (7 (Mi) and a(AMi) are the corresponding resolutions. The su- 
perscript "n" denotes the world-average mass or mass difference [Q. The sum runs over 
i = I)""; the second term in Eqn. ([!]) is not present for B° — > D + D~ candidates and 

is only present for the z = D* + term for 5° — > D* ± D T candidates. For B° — > D* + D*~ 
candidates, the average resolutions were used; for B° — > £)* ± D =F and _B° — > D + D~ de- 
cays, the resolution for each D and _D* + candidate was determined from the track co- 
variance matrices. If more than one B° _> £)(*)+£>(*)- candidate was present in a single 
event after all other selection criteria were applied, the one with the smallest Xm was 
selected. This observable is most effective for B° — > D* + D*~ since <r(AM) « 500 keV 
and 350 keV for the D°ti + final state in CLEO II and CLEO II. V, respectively. We 
require Xm < 10 f° r ^° ~~ > D* + D*~ candidates |I7|]. A typical requirement on Xm 
is < 6 for 5° -> D* ± D T (D* + -> D tt + , D° -> A-tt+tt , D" -> A+tt-tt-) and < 4 for 

5° -> D + D-(D+ -> K-7T + 7T+, D~ -> K+7T-7T-). 



C. Separation between the I? and D decay vertices 

The observable L/a(L) exploits the relatively long decay length of the D + meson (7/3cr ~ 
250 jttm) and is defined as 

L = v fl - v s ) • , (2) 

|Pd - Pel 

where vd (pd) is the reconstructed Z) candidate decay vertex (momentum). The reso- 
lution <j(L) is determined from the D candidates' covariance matrices; typically, cr(L) = 
500 (200) /iin for CLEO II (CLEO II. V). For CLEO II, only the 2-dimensional r<p informa- 
tion is precise enough to provide some discrimination so we use only the rip projection of L; 
in CLEO II. V, the SVX allows the use of the full 3-dimensional vertex information. We re- 
quire L/a(L) > for B° -> D*+D*- -> (D + ir°)(D°n-) candidates in CLEO II. V only |7 



For the CLEO II. V detector configuration, typical requirements on L/a(L) are > —0.5 
for B° -> D*±DT(D* + -> D°n+, D° -> A^Tr+Tr , D~ -> A~*~7r — 7T — ) and L/a(L) > 2.5 for 

5° -> D + D-(D + -> A^7T + 7T + , £)~ -> A+7T-7T-). 



D. Thrust and helicity angle 

For the -B — ► D^D^ and 5° — ► D + D~ modes, the observable cos #t was used to sup- 
press non-BB background. The angle between the thrust axis of the 5° candidate and 
the thrust axis of the remainder of the event is 9?. Continuum (e + e~ — > qq, q = u,c,s,d) 
backgrounds are sharply peaked towards |cos#t| — 1 and signal events are uniform in cos#t- 
The maximum allowed |cos#t| ranges from 0.50 to 0.95 for the B° — > D* ± D T channels and 
from 0.80 to 0.95 for the B° -> D+D' channels. 

The pseudoscalar — > vector, pseudoscalar decay B° — > D* ± D T produces a cos 2 #h distri- 
bution for signal and is uniform for background. The angle #h is taken between the 7r s and 
the D* + in the D* + rest frame. The minimum allowed Icos^h] for B° — ► D* ± D T candidates 
lies in the range 0.1 to 0.7, depending on the decay channel. 
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E. D + decay length 



The B° — > D + D~ mode suffers from a background that consists of a D + candidate 
where the majority of daughter candidate tracks are the result of a Z) + meson decay and 
a D~ candidate composed of a random combination of tracks. The observable L/a(L) 
(Sec. [IV Q) does not sufficiently suppress this background due to the decay length of the D + 
candidate, but a requirement on S = min(dD / 'c(^d) , dp jo (dp)), the minimum decay length 
significance of the B° daughters, where do = (vd — vg) • Pd/|pd|, reduces this background 
component. The average B° decay length is ~ 30 /im; therefore, the B° decay vertex v# 
can be accurately approximated as the e + e~ interaction point. For the K~tt + tt + , K + 7t~tt~ 
final state, we require 5* > —0.5 for the CLEO II. V configuration. 

For the B° -> D^D* and B° -> D + D~ modes, there are channels for which the back- 
ground could not be reduced to a reasonable level with any combination of selection criteria. 
Specific B° — > D* ± D T channels were discarded if the background estimated from simulation 
could not be reduced below 1/6 of the expected signal rate. Out of 84 possible channels con- 
sidered, a total of 57 and 67 B° — ► D* ± D T channels survive this criterion for the CLEO II 
and CLEO II. V detector configurations, respectively. Similarly, B° — > D + D~ channels for 
which the background could not be reduced below 1/3 or 1/7 of the expected signal rate for 
the CLEO II or CLEO II. V configuration, respectively, were rejected. These criteria select 
8 and 7 out of a total of 15 possible B° -> D+D~ channels for the CLEO II and CLEO II. V 
configurations, respectively. 



V. RESULTS AND INTERPRETATION 

The AE versus M(B) distributions of B° -> D* + D*-, B° -> D* ± D T and B° -> D + D~ 
candidates passing all selection criteria are shown in Figures [I], ^| and |3|, respectively. A 
significant signal is apparent for B° — > D* + D*~ decays; the larger backgrounds for the 
B° — ► D* ± D T and B° — > D + D~ modes are discussed below. 



A. Background estimation 

For all three modes, the background is estimated with two independent methods based 



on samples drawn largely from the data |T7|]. Method 1 uses the grand sideband (GSB) 



indicated in Figures [l], |^ and |3|. The observed number of candidates in the GSB in 
each channel is scaled to estimate the background in the signal region. The scale fac- 
tors are (7.3 ± 2.2) x 10~ 3 , (4.7 ± 1.2) xlO" 3 , (4.3 ± 0.7)xl0~ 3 and (4.0 ± 0.9)xl0~ 3 for 
the B° -> D* + D*-, B° -> D^D^CLEO II), B° -> ^^^^(CLEO II. V) and B° -> D+D- 
analyses, respectively, and are estimated from the fitted distributions in M(B) and AE. The 
excluded region of the GSB contains fully- or partially-reconstructed B — > D^*' + D^*'~X 
decays that cannot enter the signal region. The GSB regions are slightly smaller for 
the B° — > D^ + D^ analyses because they suffer from "reflection" background. "Reflec- 
tion" backgrounds arise if Cabibbo-favored B° — > £)W+£)(*)~ decays are interpreted as 
B° — ► D^ + D^ when a charged kaon from the is misidentified as a pion. This back- 
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TABLE I. Background estimates. The two background estimation methods are described in 
the text. For method 2, the combinatorial, cc and BB components of the background are listed 
separately. The uncertainties in the table are statistical only and do not include the uncertainty due 
to the background scaling factor derived from the fitted AE and M(B) distributions (Sec. VA). 



Decay 


Method 1 
Total 


Total 


Method 2 
combinatorial cc 


BB 


B° _> D*+D*~ 
B° _> D*±D=f 
B° -> D+D- 


0.384 ±0.053 
1.874 ±0.102 
0.498 ± 0.048 


0.469 ± 0.057 
1.795 ± 0.098 
0.459 ± 0.041 


0.382 ± 0.046 
1.336 ±0.062 
0.433 ± 0.039 


0.052 ± 0.034 
0.305 ± 0.078 
0.014 ±0.003 


0.035 ± 0.005 
0.064 ± 0.005 
0.013 ±0.003 



ground has AE < —50 MeV due to the kinematics of the Df decay combined with the 
difficulty in distinguishing from 7r ± for |p| > 800 MeV/c with dE/dx or TOF. 

For method 2 the contribution of each background component was estimated separately. 
The dominant contribution to the background consists of combinations of and 
in which one or both candidates is fake; that is, the daughter candidates are not the 
result of a meson decay. This combinatorial background can be estimated by forming 
explicit fake £)(*)+ candidates drawn from the D candidate mass sidebands by replacing 
Mf in Eqn. [L] with Mf + /er(Mi) or M? - fa(Mi). We use / = 6 so that classification of 
each D meson candidate as fake or standard is unique given the Xm selection criteria. The 
contribution to each channel of the combinatorial background can be derived from the two 
samples consisting of fake and standard candidates or fake Z)M + and fake D'*'" 

candidates. 

Two other background components are due to random combinations of real D^*' + 
and D"" mesons that are approximately back-to-back and arise from the processes 
e+e" -> cc -> £)W + £)W- X or e+e" -> T(4S) ^ BB ^ (D^+X)(D^-Y). The 
e + e~ — > cc — > X component was estimated from 4.6 fb _1 of e + e~ data taken 60 

MeV below the T(4S) resonance after subtraction of the combinatorial background using 
the method described above. The e+e" -> T(4S) ^ BB ^ (D^ + X)(D^~Y) component 
was estimated from samples of simulated events at least 10 times the data sample size. The 
estimated total backgrounds are listed in Table |. The estimates from the two methods for 
each channel are in good agreement and are combined channel-by-channel to produce the 
overall background estimate. 

We assess the probability for the estimated background to produce a more "signal-like" 
configuration of candidates than the observed B — > D^ + D~ signal candidates with the 
likelihood £ = Hi f [W, rii) , where the product runs over all channels selected for either the 
B° — > D* ± D T or B° — > D + D~ analysis, /(//; n) = e -M /i n /n!, b$ is the estimated background 
in the i th channel and n« is the observed number of signal candidates in the z th channel. 
We compare the distribution of C for many simulated experiments consisting solely of back- 
ground with the value of £ obtained for the signal candidates in the data. In the simulation 
of the background-only experiments, we take into account both the statistical and sys- 
tematic uncertainty in the per-channel background estimates. For the B° — > D^D^ and 
B° — > D + D~ mode, a total of 0.3% and 3.8%, respectively, of the simulated, background- 
only experiments had £ > £d&t& and, hence, are more signal-like than the observed can- 
didates. These rates are too large to claim an unambiguous observation of either the 
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B° _> D*±D=f or 5° -> D + L>-mode. For the S° -> D*+D*~ mode, fewer than 2 x 1(T 7 
background-only experiments were more signal-like than the data. 



B. Branching fraction determination 

The B° — ^ _D(*)+£)(*)~ branching fractions are determined from the likelihood 

C(B)=Y[f(^n t ) , (3) 

i 

where 

• B = B(B° -> dW+DW"), 

• fii = Si + k, 

• s 4 = 2f 00 N(BB)e i B i (DM+)B(B° -> £>M+Z>M-), 

• 6j is the reconstruction efficiency of the 2 th channel, 

• Bi(D^ + ) is the product daughter branching fractions of the z th channel and 

• N(BB) is the number of BB pairs. 

We assume f 00 /f+- = B(T(4S) -> £°5°)/£(T(4S) -> B+B~) = 1 for the results presented 
here. The evaluation of C(B) takes into account the systematic uncertainties due to the 
background estimate, efficiencies and daughter branching fractions H. The branching 

fractions and upper limits at 90% CL for the three B° decay modes are listed in Table ||. 
Since the background estimates of the two methods are combined channel-by-channel, the 
combination of the total background estimates of methods 1 and 2 (Table |) differs slightly 
from the total background estimate given in Table ||. Furthermore, the evaluation of the 
B° — ► £)(*)+£)(*) _ branching fractions with a likelihood function that takes into account 
the reconstruction efficiency, daughter branching fractions and backgrounds of each chan- 
nel (Eqn. (^)) differs from the branching fraction that would be derived from the average 
efficiency times daughter branching fraction and total backgrounds listed in Table |I[ 

While only the B° — > D* + D*~ results provide unambiguous evidence of the Cabibbo- 
suppressed h —>■ ccd decay, the expectations based on the corresponding Cabibbo-favored 
decays are consistent with the upper limits of the other two modes. The results presented 
here indicate that there may be potential difficulties in the measurement of sin 2(3 using 
B° — > Z}(*)+_D(*) _ decays. The yields are appreciably lower than that of B° — > J/ipKg for the 
same integrated luminosity, and background levels are higher, especially for B° — ► /}* ± D =F 
and B° — > D + D~ . Measurement of sin 2(3 via the proper-time dependence of B° — > D^ + D~ 
decays performed at asymmetric e + e~ colliders or at hadron colliders may be able to exploit 
the B° decay length to reduce backgrounds. In contrast, the B° — > D* + D*~ results show 
that this mode, while also having a yield substantially lower than that of B° — > J/ipKg, has 
very low backgrounds and should provide an independent measure of sin 2(3. The suppression 
of background for B° — > D* + D*~ is achieved largely through the observable Xm(Scc. IV B ) 



that relies on accurate reconstruction of the trajectory of the charged slow pion from the D* 
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decay. Inability to reconstruct efficiently the 7r+ can substantially degrade a potential sin 2/3 
measurement. For example, for the results presented here, the reconstruction efficiency 
of the 7r+ from D* + — > P°7r+ for the CLEO II. V configuration is (65 ±6)% of that for the 
CLEO II configuration because the track-finding algorithm was optimized only for the latter 



configuration fll? 



TABLE II. The number of observed candidates, estimated total backgrounds, efficiencies, 
measured branching fractions and branching fraction upper limits at 90% CL for the three 
B° — > Z)(*) + £)(*)"" modes. For the branching fractions and background, the first error is the 
statistical uncertainty and the second is the systematic uncertainty. (eB) is the product of the 
reconstruction efficiencies and the daughter branching fractions summed over all channels; 
the uncertainty includes both the statistical uncertainty in the estimation of e, from simulation as 
well as the uncertainties in the daughter branching fractions ||. 



Decay 
mode 




Candidates 


Total 
background 


(eB) 
(xlO" 4 ) 


Branching 
fraction (xlO -4 ) 


90%CL Upper 
limit (xlO" 4 ) 


B ^ 
B° — > 
B ^ 


D*+D*~ 

D*±DT 

D+D- 


8 
6 
2 


0.42 ±0.04 ±0.13 
1.68 ±0.07 ±0.24 
0.46 ±0.03 ±0.10 


8.2 ±2.9 
11.0 ± 1.7 
5.4 ± 1.0 


9.9j£jj ± 1.2 
2.1± 2 ; 4 .±0.5 

o n +3. 3+0.8 
- u -2. 1-0.6 


6.3 
9.4 



C. B° — » D* + D* transversity analysis 

A measurement of sin 2/3 from B° — > D* + D*~ decays requires an angular analysis to 
disentangle the CP-odd and CP-even components of the decay. In the transversity basis |TU| 



the fraction of the CP-even component (A) of the decay B° — > D* + D* can be determined 
from the cos 9 tr distribution, 

' ' /J ^sin 2 ^.±^(l-A)cos 2 ^ r , (4) 



Tdcos6 tr 4 ' 2 



where T = T(B° D* + D*~) + T(B° -> D* + D*-) and 6 tr is the angle between the vr s from 
the D* + and the normal to the plane of the D*~ decay in the D* + rest frame as shown in 
Fig-! 

We perform an unbinned, maximum likelihood fit to extract A from the cos 6^ distri- 
bution of the eight B° — > D* + D*~ candidates, taking into account the background shape 
and the cos 6^ resolution and acceptance. The background shape, estimated from GSB 
candidates, is consistent with being uniform as a function of cos 6^. The resolution of 
<r(cos#t r ) = 0.1 is determined from simulated events in the observed decay channels. The 
acceptance varies as a function of cos 9 tr due to the drop in efficiency at low momentum for 
the charged 7r s . For 7r+ emitted perpendicular (parallel) to the D* + direction, cos 9 tr tends 
towards ±1 (0). Thus a loss of efficiency for low momentum 7r+ results in a reduction of 
acceptance at cosd tr near zero. This effect is inconsequential for the D* + — ► D + tt° candi- 
dates because the 7r° efficiency does not vary appreciably. The acceptance is modeled as 
oc 1 + a cos 2 9 tr , where a — 0.17 ± 0.17 is determined from simulated B° — > D* + D*~ decays 
and the uncertainty represents a conservative estimate of the range of a. 
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The observed cos 9t r distribution of the B° — > D* + D*~ candidates is shown in Fig. || with 
the fit result superimposed. Figure || shows the dependence of L(A) = —2hi(C(A)/C(A)), 
assuming a = 0.34 where A is the value of A that maximizes C(A). The conventional 
evaluation of confidence levels from L(A) is confounded because the statistical resolution on 
A is comparable to the bounds on A of [0, 1]. To determine confidence levels, we evaluate 
L(A) as a function of the input value of A using 10000 simulated experiments at each value 
of A n p Ut = 0.0, 0.1, 0.2, . . . , 1.0. Each simulated experiment is analyzed as the data and the 
distribution of dN/dL(A inpnt ) is determined (N is the number of simulated experiments). 
At each value of A input , we then determine the 95% CL value, L 95 , as 

L 95 dN I r°° , dN 

dL — dL — = 0.95 . 5 
o dL 1 Jo dL 

In Fig. |6] we show the curves resulting from this procedure at the 68.3, 90, 95 and 99% 
CL for a = 0.34. The confidence level curves have a concave shape because the dN/dL 
distributions peak more sharply at A input near and 1 due to the bounds on A. We perform 
this procedure for the central and extreme values of the acceptance, a = 0.00,0.17,0.34, 
for both the simulation and the data to take into account the acceptance uncertainty. We 
conservatively use the regions excluded by all three values of a to set limits. We exclude 
values of A < 0.11 at 90% CL, but cannot exclude A = at 99% CL. Combining the limits 
0.15 < A < 0.90 at 68.3% CL for the three values of a with the most likely value of A for 
a = 0.17 and taking into account the uncertainties in the level and shape of the background, 
we find A = 0.49to;34 ± 0.02. Our results are consistent with expectations that A « 0.95 
although the statistical precision is poor. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the decays B° -> D* + D*~, B° -> D*^D T and B° -> D+DAn 9.7 x 10 6 
T(4S) -> BB decays. We determine B(B° -> D* + D*-) = (9.9±|;1 [stat.] ± 1.2 [syst.]) x 10~ 4 
and limit B(B° -> D^D^) < 6.3xl0~ 4 and B(B° -> D+D~) < 9.4xl0~ 4 at 90% CL. These 
results, while consistent with expectations, show that substantially higher luminosities will 
be needed to make a measurement of sin 2(3 using B° — > decays that approaches 

the statistical precision of a sin 2(3 measurement using B° — > J/ipK®. Asymmetry mea- 
surements of lesser precision with B° — > Z}(*)+Z}(*) _ decays may, however, be adequate for 
resolving ambiguities in the determination of (3. We have performed the first transversity 
analysis for B° — > D* + D*~ and exclude values of the CP-even component of the decay less 
than 0.11 at 90% CL. 
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FIG. 1. (a) The AE vs. M(B) distribution for B° -> D* + D*~ candidates for the data taken 
at the T(4S) resonance. The small rectangle delineates the signal region and the region outside 
the dashed line is the GSB. (b) The M(B) distribution with the requirement \AE\ < 20 MeV. 



15 



0.4 



3080100-004 



0.2- . 



~ 



-0.2 



>-0.4 



(a) 



• •• 



*. ... 



. • y 



(b) 



5.20 



5.22 



5.24 5.26 
M(B) (GeV) 



5.28 



FIG. 2. (a) The AE vs. M(B) distribution for B° -» D* ± D =F candidates for the data taken 
at the T(4S) resonance. The small rectangle delineates the signal region and the region outside 
the dashed line is the GSB. (b) The M{B) distribution with the requirement \AE\ < 16 MeV. 
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FIG. 3. (a) The AE vs. M(B) distribution for B° -> D + D~ candidates for the data taken at 
the T(4S) resonance. The small rectangle delineates the signal region and the region outside the 
dashed line is the GSB. (b) The M{B) distribution with the requirement \AE\ < 16 MeV. 
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FIG. 5. The fitted cos 6f r distribution of the eight B° — > D* + D*~ candidates from the signal 
region. The filled histogram represents the data, the solid line represents the best fit result and 
the dashed line represents the background component. The fit takes into account the acceptance 
and resolution in cos 6f r as described in the text. 
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A = Fraction of CP-even Component 

FIG. 6. L data (A) = -2\u{C{A) / C{A)) for the data (solid curve) compared to L n (Eqn. ©) 
for n = 68.3%, 90%, 95% and 99% (broken lines) that correspond to the confidence levels at n% 
for a = 0.34. See text for details. 
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